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y—{ Abstract. In this paper wc study the geometric discrepancy of explicit constructions of 

uniformly distributed points on the two-dimensional unit sphere. We show that the spherical 
cap discrepancy of random point sets, of spherical digital nets and of spherical Fibonacci 
lattices converges with order iV~^/^. Such point sets are therefore useful for numerical 
^ ^ integration and other computational simulations. The proof uses an area-preserving Lambert 

^ map. A detailed analysis of the level curves and sets of the pre-images of spherical caps under 

this map is given. 

'"7 1. Introduction 

^ Let = {z G : \\z\\ = 1} be the unit sphere in the Eudidean space provided with 

the norm || ■ || induced by the usual inner product x ■ y. On this sphere we consider the 



Lebesgue surface area measure a normalised to a probability measure (/g2 do" = 1) 



^ This paper is concerned with uniformly distributed sequences of points on Informally 

^ G ^ speaking, a sequence of points is called uniformly distributed if every reasonably defined 

(clopcn) y4 C §2 gets a fair share of points as their number N grows. Given a triangular 
^ scheme {zi^n, . . . , zn,n}, N > 1, of points on in such a case one has 

^ lim card({j : Zj, n £ A}) _ 

^ (where card denotes the cardinality of the set) or, equivalently (defined in terms of numerical 

integration) 

O 1 r 

^ (lb) lim -r;/^f{zjN) = / /d<7 for every / continuous on 

• • 7 = 1 

> 

• ^ The degree of uniformity is quantified by the so called spherical cap discrepancy. 

^ A spherical cap C = C{w,t) centred at G §^ with height t E [—1, 1] is given by the set 

<^ C{w, t) = {yeS'^ -.w-yyt}. 

(We assume that spherical caps are open subsets of S^.) The boundary of C{w,t) then is 

dC{w,t) = {yeS'^:w-y = t}. 
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Let Z]\f = {zq, . . . , ^AT-i} C §2 be an A^-point set on the sphere S^. The local discrepancy 
with respect to a spherical cap C measures the difference between the proportion of points 
in C (the empirical measure of C) and the normalised surface area of C. The spherical cap 
discrepancy is then the supremum of the local discrepancy over all spherical caps, as stated 
in the following definition. 

Definition 1. The spherical cap discrepancy of an A^-point set = {zq, . . . , z^-i} C is 



D{Zn) 



sup 

K;eS2 ■ 



sup 

■ l<t<l 



a{C{z,t)) 



n=0 



If the point set Z^r is well distributed, then this discrepancy is small. In fact, a sequence 
of A^-point systems (Z7v)Ar>i satisfying 



(2) 



lim D{Zn) 



0, 



is called asymptotically uniformly distributed. Using, for example, the classical Erdos-Turdn 
type inequality (cf. Grabner f^S], also cf. Li and Vaaler [38]) or LeVeque type inequalities 
(Narcowich, Sun, Ward, and Wu |41j ) and the fact that the set of polynomials is dense in the 
set of continuous functions, one can show that ^ is equivalent with (lb). 

It is known from [6J that there are constants c,C > 0, independent of N, such that a 
low- discrepancy scheme {Z'^}n>2 satisfies 



(3) 



The lower bound holds for all A^-point sets Zjy on and there always exists an A'^-point set 
Zjv C §2 such that the upper bound holds. The proof of the upper bound is probabilistic 
in nature and is thus non-constructive. To our best knowledge explicit constructions of low- 
discrepancy schemes are not known. (In this paper we restrict ourselves to the sphere S^, 
though some of the results are known for spheres of dimension d >2.) 

An explicit construction of points Z^- with small spherical cap discrepancy has been given 
in \39\ Ho] . For instance, in [39] it was shown that 

(4) D{Zn) <C{logN)^^^N-^/^. 

The numerical experiments in [39] indicate a convergence rate ofO{N~^/'^). 
In this paper we give explicit constructions of point sets Zf^ for which we have 

D{Zn) < 44\/2A^-i/2^ 

Our numerical results indicate a convergence rate of 0{{log NyN-^/^) for some 1/2 < c < 1, 
see Tables [1} [2] below. 

The spherical cap 1^2- discrepancy 

1/2 



Di.,{Zn) 




N-l 



a{C{z,t)) 



n=0 



da{w) dt 



which averages the local discrepancy for a spherical cap over all caps, provides a lower bound 
for the spherical cap discrepancy. It is closely related to the sum of distances and its continu- 
ous counterpart the distance integral by means of Stolarsky 's invariance principle [50] for the 
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Euclidean distance and the 2-sphere, 

N N „ „ . 

]^EEl-"j--^'^l+4K,(^i,---,^7v)]'= / / \z-w\da{z)da{w)=:V^,{S') = -. 
j=i k=i ''^^ -^^^ 

This gives a simple way of computing the spherical cap L2-discrepancy of point sets on S^. 
In [H] it is shown that the spherical cap L2-discrepancy of Z]\f can be interpreted as the 
worst-case error of an equal weight numerical integration rule with node set Zj^ for functions 
in the unit ball of a certain Sobolev space over S^. It is shown in [13j that on average (that is 
for randomly chosen points independently identically uniformly distributed over the sphere), 
the expected squared worst-case error is of the form (4/3) A^^^. Thus the expected value of 
the squared spherical cap discrepancy satisfies 

(5) 8E [D{ZN)f > 4E [Di^,{ZN)f = ^ N-\ 

We study the expected value and the typical asymptotic order of the spherical cap discrepancy 
of random point sets in detail in Section [4} Among other results, we show that the there is 
also a constant C > such that E [D{Zn)] < CN^^/'^. 



Point configurations maximising the sum of distances, by Stolarsky's invariance princi- 
ple, have low spherical cap L2-discrepancy. It is known from [6] that low spherical cap 
L2-discrepancy point sets satisfy relations similar to ([s]) except for the logarithmic term in- 
troduced by the probabilistic approach. The upper bound for the spherical cap discrepancy 
of maximum sum of distances points Z"^ obtained in [H] is much weaker but still better than 
For some positive constant c > 0, not depending on A'^, 

D{Z*^)< cN'^^^. 

For point configurations Z'^ emulating electrons restricted to move on in the most stable 
equilibrium, that is minimising their Coulomb potential energy essentially given by 

N N 



EE 



1 



i=lfc=l'--^ ^^=1 

one can show the bound 

D{Z*j^) < CN^^/^ log N. 

The estimate D{Z*j^) = 0{N~^/'^) was conjectured by Korevaar [35] and later proved (up to 
the logarithmic factor) by Gotz [21]. When allowing so-called K -regular test sets [^j introduced 
by Sjogren [3^], the estimate above is sharp in the following sense: The upper bound holds 
for any ET-regular test set, whereas there are some numbers Kq and c such that to any N 
points zi, . . . , zjy S there is a ^Co-regular test set B with pi's Corollary 2] 

1 ^ 

n=l 

(The lower bound also applies to the explicit constructions given in this paper.) Bounds 
for the spherical cap discrepancy of so-called minimal Riesz energy configurations (for the 



'Roughly speaking, cr-measurable subsets B of whose (5-neighbour hoods dsB relative to satisfy 
a{dsB) < KS. For example, spherical caps are -fTi-regular for some fixed Ki. 
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concept of Riesz energy see, for example Saff and Kuijlaars [35] and Hardin and Saff |28j ) 
can be found in (for the logarithmic energy), Damelin and Grabner \17[ [T8] (the first 
hyper-singular case), and |41j (sums of generalised distances). Wagner [53] estimates the 
spherical cap discrepancy in terms of the Riesz energy. It should be mentioned that there are 
very few known explicit constructions of point configurations with optimal Riesz energy. In 
general, one has to rely on numerical optimisation to generate such point sets. The underlying 
(constrained) optimisation problem is highly non-linear. Moreover, numerical results indicate 
that the number of local minima increases exponentially with the number of points. (For the 
computational complexity see, for example, Bendito et al. [9].) 

Spherical n-designs introduced by Delsarte, Goethals and Seidel in the landmark paper [19j 
are node sets for equal weight numerical integration rules such that all spherical polynomials 
of degree < n are integrated exactly. Grabner and Tichy [26] give the following upper bound 
of the spherical cap discrepancy of a spherical n-design with N{n) points 



which immediately follows from the aforementioned Erdds-Turdn type inequality. (See also 
Andrievskii, Blatt, and Gotz |4] for a similar form for X-regular test sets.) 

A spherical n-design is the solution of a system of polynomial equations (one for every 
spherical harmonic of the real orthonormal basis of the space of spherical polynomials of 
degree < n). Hence, a natural lower bound for the number of points of a spherical n- 
design is given by the dimension of the involved polynomial space; that is, one needs at 
least > points. The famous conjecture that C points (for some universal C > 0) 

are sufficient for a spherical n-design seems to have been settled by Bondarenko, Radchenko 
and Viazovska [ITO]. The proposed proof is non-constructive. Hardin and Sloane [29] propose 
a construction of so-called putative spherical n-designs with (1/2) n^ -|- o(n^) points. The 
variational characterisation of spherical designs introduced in [IH] (also cf. |26j ) leads to a 
minimisation problem for a certain energy functional (changing with n) whose minimiser is a 
spherical n-design if and only if the functional becomes zero. Numerical results also suggest a 
coefficient 1/2. When allowing more points, N{n) = (n-|-l)'^, interval-based methods yield, in 
principle, the existence of a spherical n-design near so-called extremal (maximum determinant 
points, cf. [?7j)- Due to the computational cost this approach was carried out only for n < 20. 
Very recently, Chen, Frommer, and Lang [16j devised a computational algorithm based on 
interval arithmetic that, upon successful completion, verifies the existence of a spherical n- 
design with (n-f-l)^ points and provides narrow interval enclosures which are known to contain 
these nodes with mathematical certainty. The spherical cap discrepancy of all such obtained 
spherical n-design with 0{'n?) points can then be bounded by C N~^/^ by ([6]). For the sake 
of completeness it should be mentioned that the tensor product rules used by Korevaar and 
Meyers [3^ to prove the existence of spherical n-designs of A^(n) = 0{n^) points give rise to 
A^(n)-point configurations whose spherical cap discrepancy can be bounded by C" [A^(n)]~^/^ 



From [H] it follows that the spherical cap discrepancy of a point set Zj^ = {zq, . . . , ^tv-i} ^ 
yields an upper bound on the integration error in certain Sobolev spaces of functions defined 
on using a quadrature rule Qnif) = Yln=o fi^n)- Thus, our results here provide an 
explicit mean of finding quadrature points for numerical integration of functions defined on 
S^. Our result here improves the bound on the integration error in [12j by a factor of \/\ogN. 

The construction of the points on is obtained by mapping low-discrepancy points on 
[0, 1]2 to §2 using an equal area transformation $ : [0,1]^ — )• S^. The same approach has 



(6) 




by (g. 
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previously been used in |12| and [27j, in both cases in the context of numerical integration. 
The low-discrepancy points in [0, 1]^ are obtained from digital nets and Fibonacci lattices, 
see [201 m]. These point sets are well-distributed with respect to rectangles anchored at the 
origin (0,0). However, the set 

^-\C{w, t)) = {xG [0, 1]2 : ^{x) £ C{w, t)} 

is, in general, not a rectangle. In fact, it is not even a convex set, although the boundary 

^-\dC{w, t)) = {xe [0, 1]^ : ^{x) G dC{w, t)} 

is a continuous curve. 

Hence, in order to prove bounds on the spherical cap discrepancy of digital nets and Fi- 
bonacci lattices lifted to the sphere using $ (we call those point sets spherical digital nets 
and spherical Fibonacci lattices), we need to prove bounds on a general notion of discrepancy 
in [0, 1]^. To this end we study discrepancy in [0, 1]^ with respect to convex sets, the corre- 
sponding discrepancy is known as isotropic discrepancy [8]. We show that digital nets and 
Fibonacci lattices have isotropic discrepancy of order 0{N~^/'^). Using these result and some 
properties of the function we can show that spherical digital nets and spherical Fibonacci 
lattices have spherical cap discrepancy at most CN~^/'^ for an explicitly given constant C, 
see Corollary [16] and [Tsj Note that the best possible rate of convergence of the isotropic 
discrepancy is A^~^/^(log A^)'^ for some < c < 4, see [7] and [8^, p. 107]. Hence the approach 
via the isotropic discrepancy cannot give the optimal rate of convergence for the spherical 
cap discrepancy. 

In the following we define the equal area Lambert map $ and show some of its properties. 

2. The equal-area Lambert transform and some properties 
The points on the sphere are obtained by using the Lambert cylindrical equal-area projection 
(7) $(a, r) = {2^Jt - cos(27ra), 2\/r - sin(27ra), 1 - 2t) , a,rG [0,1]. 

The area-preserving Lambert map can be illustrated in the following way. The unit square 
[0, 1]^ is linearly stretched to the rectangle [0, 27r] x [—1, 1], rolled into a cylinder of radius 1 
and height 2 and fitted around the unit sphere such that the polar axis is the main z-axis. This 
way a point (a,T) in [0, 1]^ is mapped to a point on the cylinder which is radially projected 
along a ray orthogonal to the polar axis onto the sphere giving the point $(a,T). 

Axis-parallel rectangles in the unit square are mapped to spherical "rectangles" of equal 
area, see Figure [ij 

The pre-images of a spherical cap centred at w with height t under the Lambert map is 
the set 

B{w, t) = ^-^{C{w, t)) = {(a, r) G [0, 1) x [0, 1] : *(a, r) G C{w, t)] 

and the pre-image of the boundary of this spherical cap is 

dB{w, t) = i^'^{dC{w, t)) = {{a, r) G [0, 1) x [0, 1] : *(a, r) G dC{w, t)] . 

The sets B{w, t) are not convex, in general. Thus, we consider a more general class of sets 
which we call pseudo-convex. A definition is given in the following. 

Definition 2. Let A be an open subset of [0, 1]^ such that there exists a collection of p convex 
subsets Ai, Ap of [0, 1]^ with the following properties: (a) Aj n Ak is empty for j / k; 
(b) A C AiD ■ ■ ■ U Ap] (c) either Aj is a convex part of A [Aj C A) or the complement of 
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A with respect to Aj, A'- = Aj \ A, is convex. Then A is cahed a pseudo-convex set and 
Ai, . . . , Ap is an admissible convex covering for A with p parts (with q convex parts of A). 

Lemma 3. For every w € and all —1 < t < 1 the pre-image B{w, t) of the spherical cap 
C{w, t) centred at w with height t under the Lambert map is pseudo-convex with an admissible 
convex covering with at most 7 parts. More precisely, taking into account the number of convex 
parts of the pre-image, among the convex coverings with p parts and q of which are convex 
the worst case has p = 7 and q = 3 which implies the constant 2p — q = 11. 

The proof of Lemma |3] in Section [7] gives details how to construct admissible coverings. 

3. Isotropic- and spherical cap discrepancy 

We introduce the isotropic discrepancy of a point set and a sequence as follows. Let A be 
the Lebesgue area measure in the unit square. 

Definition 4. The isotropic discrepancy Jn of an A^-point set Pn = {a;o, . . . , a^Tv-i} in [0, 1]^ 
is defined as 

N-l 



N) 



sup 



1 

N ^ 

n=0 



lA{Xn) - XiA) 



where A is the family of all convex subsets of [0, 1]^. 

For an infinite sequence xq,xi,... S [0,1]^ the isotropic discrepancy is defined as the 
isotropic discrepancy of the initial N points of the sequence. 

Lemma 5. Let A be a pseudo-convex subset of [0, 1]^ with an admissible convex covering of 
p parts with q convex parts of A. Then for any N -point set P/v = {xq, . . . , a^TV-i} ^ [0, 1]^ 

N-l 

- J2 UK) - A(^) 



n=0 



< {2p-q)JN{PN). 
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Proof. Let Ai, . . . , Ap be an admissible convex covering of A with p parts. Without loss of 
generality let , . . . , be the convex parts of A and Aqj^i , . . . , yip those for which A'- = Aj\A 
(g + 1 < J < p) is convex. Clearly 

Thus 

1 

- UK) - KA) 



n=0 



Af-1 



MA, 



j=l n=0 
j=l n=0 



N-l 



j=q+l n=0 
P ^ N-l 



j=q+l 



n=0 



j=g+l n=0 

In the last line all sets are convex and we can use the isotropic discrepancy in the estimation 

Af-l 

< [q + {p - q) + (p - q)] Jn{xo, xn-i). 



-YlA{Xn)-X{A) 



n=0 



□ 



Theorem 6. Let Pn = {xq, . . . , a^Tv-i} ^ [0, 1]^ and let = {^{xq), . . . , ^{x^^i)} C S"^. 



Then 



D{Zn) < IIJn{P. 



N) 



Proof. Let w € and —l<t<l. A point $(a;„) G C{w,t) if and only if a;„ € B{w,t). 
Thus 

N-l N-l 
n=0 n=0 

Further, since the transformation $ preserves areas, we have 



a{C{w,t)) = X{B{w,t)). 



Hence 



N-l 



N-l 



^Y^C:iw,t){^{Xn))-(j{C{w,t)) = ^Y,^B{w,t){Xn) - X{B{w,t)) 
n=0 n=0 

The pre-images are pseudo-convex in the sense of Definition[2]by Lemma[3j Applying Lemma[5] 
with the constant 2p — q = 11 from Lemma [2] we arrive at the result. 



□ 



We have now reduced the problem of proving bounds on the spherical cap discrepancy to 
prove bounds on the isotropic discrepancy of points in the square [0, 1]^. We will study this 
problem in Section [5} 
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4. Spherical cap discrepancy of random points sets 

Let {M,Ai) be a measurable space, and let P be a probability on it. Let further X„, 
n > 0, denote a sequence of independent, identically distributed (i.i.d.) random variables on 
a probability space (fi,^, P) with values in M, and let C C denote a class of subsets of 
M. To avoid measurability problems we will assume throughout the rest of this section that 
the class C is countable. Let A C M be an arbitrary set. Then C is said to shatter A if to 
every possible subset B of A there exists a set C G C such that 

CnA = B. 

For k > 1 the k-th shattering coefficient Sc{k) of C is defined as 

Sc{k):= max card{{xi, . . . , x^} CiC : C £ C}. 

Xi,...,Xfe£A/ 

The Vapnik-Cervonenkis dimension (VC-dimension) of C is defined as 

v{C):=mm{k : Sc < 2^}. 

k 

(Here we use the convention that the minimum of the empty set is oo.) A class C with finite 
VC-dimension is called a Vapnik- Cervonenkis class ( VC class) . The theory of VC classes is of 
extraordinary importance in the theory of empirical processes indexed by classes of functions. 
For example, a class C is uniformly Glivenko-Cantelli if and only if it is a VC class, see 
[53]. We will use the following theorem, which is a combination of results of Talagrand |5H 
Theorem 6.6] and Haussler [30i, Corollary 1], and has already been used by Heinrich et al. 
|31j in the context of probabilistic discrepancy theory. 



Theorem 7 (see |3H Theorem 2]). There exists a positive number K such that for each VC 
class C and each probability P and sequence X„, n > 0, as above the following holds: For all 
s > K^Jv{C) we have 



sup 

cec 



- lc(^n) - P{C) 

n=0 



1 fKs' 



> ^ > < - ^ e 



N s \v{C) 



v{C) 



'2 



,2 



In our setting we will have M = 'EP' , M will denote the sigma- field generated by the class of 
spherical caps, P will stand for the normalised Lebesgue surface area measure a, and C will 
denote the class of all spherical caps for which the centre w is a vector of rational numbers 
and the height t is also a rational number (this restriction is necessary to assure that the class 
C is countable; of course the spherical cap discrepancy with respect to this class is the same 
as the discrepancy with respect to the class of all spherical caps). In the sequel we assume 
that the i.i.d. random variables X„, n > 0, are uniformly distributed on S^. We will write 
Zjv = Zn{u}) for the (random) point set {Xq, . . . , Xjy^i} = {Xo{u}), . . . , XAr_i(u;)}. 

The following proposition asserts that the class C is a VC class (the proof of this and the 
subsequent results of this section can be found in Section [7]) . 

Proposition 8. The class C has VC dimension 5. 

Using Theorem [7] and Proposition [8] we can prove the following results: 

Theorem 9. There exist constants Ci,C2 such that for N >1 

CiVn < E [D{Zn)] < CaViV. 
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Remark. The existence of such a constant Ci for the lower bound fohows directly from ([5]); 
we can choose Ci = 6^^/^. 

Theorem 10. For any e > there exist positive constants C3{e),C4{e) such that for suffi- 
ciently large N 

F |C73 < y/ND{ZN) < C4} > 1 - e. 



Theorem 10 shows that the typical discrepancy of a random set of points is of order 
A^~^/^. However, actually much more is true, since by classical results any VC class C on a 
measurable space (M, A^) is a so-called Donsker class, which essentially means that for every 
probability measure P and every sequence Vn, n > 0, of i.i.d. random variables having law P 

the empirical process indexed by sets 

1 

-yic{Vn)-p{c) 



aN{C) = VN 



N ^ 

n=0 



C eC 



converges weakly to a centered, bounded Gaussian process B(C), which has covariance struc- 
ture 

EB{Ci)B{C2) = P{Ci n C2) - P{Ci)P{C2), Ci, C2 G c. 

This weak convergence could, for example, be used to prove the existence of a limit distribution 
of VnD{Zn) as — )• 00; however, to keep this presentation short and self-contained we will 
not pursue this method any further, and refer the interested reader to [21 [211 1221 [S2] and the 
references therein. 

Remark. The upper bounds in Theorem [9] and [To] follow from Theorem [7[ However, since 
no concrete value for the constant K in Theorem [7] is known, the value of the constants C2 
and C4 in Theorem [9] and Theorem 10, respectively, is also unknown. It is possible that the 



decomposition technique from [T] can be used to achieve a version of Theorem [9] and 10 with 
explicitly known constants in the upper bound. 

Finally, the following theorem describes the asymptotic order of a typical infinite sequence 
of random points. 

Theorem 11. We have 

D{Zn) = O ^JJ^ — 1 as N ^ 00, almost surely. 



Theorem [TT] is a so-called hounded law of the iterated logarithm, and follows easily from 
Theorem |6] and Philipp's law of the iterated logarithm (LIL) for the isotropic discrepancy of 
random point sets in the plane. More precisely, Philipp [H] proved that for a sequence of 
i.i.d. uniformly distributed random variables 1^, n > 0, on the unit square (writing P/v for 
the (random) point set {Yq, . . . , YAf-i}) the law of the iterated logarithm 

NJn{Pn) 1 
limsup , = - a.s. 

iv^oo V2iVloglogiV 2 

holds. Together with Theorem [g] this implies for = {*(lo), • • • , *(^Af-i)} ^ 

ND{Zn) ^ 11 
iim sup , < — a.s., 

Tv^oo V2A^loglogiV - 2 
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which proves Theorem 11 (it is necessary to observe that the image of a sequence of i.i.d. 
uniformly distributed random variables on the unit square under the area-preserving Lambert 
map is a sequence of i.i.d. uniformly distributed random variables on the sphere). It is easy to 



see that Theorem 11 is optimal, except for the value of the implied constant. More precisely, 
let C* denote a fixed spherical cap with area 2ii (which means that C* is a hemisphere, and 
has normalised surface area measure o"(C*) = 1/2). Then clearly the random variables 

lc-(*(r„))-a(C*), n>0, 

have expected value and variance 1/4. Thus by the classical law of the iterated logarithm 
for sequences of i.i.d. random variables 



lim sup 

Af-s>oo 



^En=0 lc*(*(>^n))-l/2 



\/27Vlog \ogN 



a.s. 



and since 



we finally arrive at 



D{Zn) > 



lim sup 

N- 



N ^ 

n=0 



(*(y„))-a(C7*) 



> 



a.s. 



V2iVlog logN - 2 

which proves the optimality of Theorem [Tl] We remark that it should also be possible to 
prove Theorem 1 1 1 1 without using Theorem[6]and Philipp's LIL for the isotropic discrepancy, by 
deducing it directly from the bounded LIL for empirical processes on VC classes of Alexander 



and Talagrand [3]. We conjecture that Theorem 11 can be improved to 

ND{Zn) _ _ 1 
^ ~ 2 



lim sup 



N^oo V2iVloglogAf 
but this seems to be very difficult to prove. 



a.s. 



5. Point sets with small isotropic discrepancy 

In this section we investigate the isotropic discrepancy of (0, m, 2)-nets and Fibonacci 
lattices. In particular we show that the isotropic discrepancy of those point sets converges with 
order O^N'^^"^). Note that the best possible rate of convergence of the isotropic discrepancy 
is 

Ar-2/3(iogAr)c for some < c < 4, see [7] and j8l p. 107]. Whether (0, m, 2)-nets and/or 
Fibonacci lattices achieve the optimal rate of convergence for the isotropic discrepancy is an 
open question. 

5.1. Nets and sequences. We give the definition of (0, m, 2)-nets in base b in the following. 

Definition 12. Let b > 2 and m > 1 be integers. A point set P;,™ C [0, 1)^ consisting of 6"* 
points is called a (0, m, 2)-net in base b, if for all nonnegative integers di,d2 with di + d2 = m, 
each of the elementary intervals 



n 

i=l 



+ 1 



< Oj < b'^^ {ai an integer). 



contains exactly 1 point of P^m. 
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It is also possible to construct nested (0, m, 2)-nets, thereby obtaining an infinite sequence 
of points. 

Definition 13. Let 6 > 2 be an integer. A sequence xq,xi,... G [0,1)^ is called a (0,2)- 
sequence in base b, if for all m > and for all A; > 0, the point set Xkb^, Xkb"^+i, ■ ■ ■ , ^{k+i)b"'-i 
is a (0, m, 2)-net in base b. 

Explicit constructions of (0, m, 2)-nets and (0, 2)-sequences are due to Sobol' [49j and 
Faure [23], see also [20l Chapter 8]. 

The following is a special case of an unpublished result due to Gerhard Larcher. For 
completeness we include a proof here. 

Theorem 14. For the isotropic discrepancy Jn of a (0, m, 2) -net P/v in base b [N = b"^) we 
have 

JAr(P^)< 4^/26- LW2J < 



iV 

Proof. Let P/v = {a^o, • • • , x^m^i}. Let k = [m/2j and consider a subcube W of [0, 1)* of the 
form 



W 



Cl Cl + 1 



C2 C2 + 1 



with ^ < Ci < b^ {ci an integer) for i = 1,2. The cube W has volume b "^^ and is the union 
of b"^~'^^ elementary intervals of order m. Indeed, 

^m — 2fc 



W 



U 

v=0 



Cl Cl + 1 



C2 



+ 



C2 V + 1 



So W contains exactly b"^~'^^ points of the net. The diagonal of W has length y^/i'^. 

Let now A be an arbitrary convex subset of [0, 1]^. Let W° denote the union of cubes W 
fully contained in A and let W denote the union of cubes W having non-empty intersection 
with A or its boundary. The sets W and W° are fair with respect to the net, that is, 

N-l N-1 

-^%K) = A(iy) and -Y^lw.{xr.) = X{wn. 



n=0 



n=0 



We have 



Af-l 



Af-1 



1 UK) - A(^) <^J2 ^wM - KW) + \{W \ A) = X{W \ A) 



n=0 



n=0 



and 



N-l 



N-l 



1 ^ UiXn) - HA) > ^ E - ^(^°) - \ = -^(^ \ ^°)- 

n=0 n=0 

Since the set A is convex, the length of the boundary of A is at most the circumference of 
the unit square, which is 4. Further we have 

W\A'Z{xe[0,lf\A: \\x - y\\ < Vlb^'' for some y e A} 

and therefore 

X{W \A)< X{{x G [0, lf\A: \\x - y\\ < ^26"'= for some y e A}) < 4^/26"'=, 
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where the last inequaUty follows from the fact that the outer boundary of the enclosing set 
has length at most 4 (which is the circumference of the square [0, 1]^). Moreover 

A\W° <^ {x A : \\x - y\\ < Vlb"^ for some y G [0, 1]^ \ A} 

and therefore 

\ W°) < X{{x €A:\\x- y\\ < ^26"^ for some y G [0, 1]^ \ A}) < 4^/26~^ 

since, by the convexity of A, the boundary of A has length at most 4 (which is the circum- 
ference of the square [0, 1]^). 
Thus we obtain 

N-1 

- ^ U(a;„) - A(^) <4V2r^ 

n=0 

and hence the result follows. □ 

Note that the above result only applies when the number of points is of the form N = 
(notice that choosing m = 1 only yields a trivial result, hence one usually chooses a small 
base b and a 'large' value of m). In the following we give an extension where the number of 
points can take on arbitrary positive integers. 

Theorem 15. For the isotropic discrepancy J]\j of the first N points P/v = {xq, . . . , a^TV-i} 
of a (0, 2) -sequence in base b we have 



JN{PN)<4V2{b' + b^/^)^ 



N 

Proof Let iV G N have base b expansion N = No+Nib-\ hNmb"^. Let Pn = {xq, . . . ,xn-i} 

denote the first N points of a (0, 2)-sequence in base b. For < k < m with > and 
< £ < iVfc the point set 

is a (0, k, 2)-net in base b by Definition [l3] Thus P/v is a disjoint union of such (0, k, 2)-nets 



^N= \J U Qk, 

0<k<mO<e<Nk 
Affc>0 



We have the following triangle inequality for the isotropic discrepancy (which is an analogue 
to the triangle inequality for the star-discrepancy [SI P- 115, Theorem 2.6]) 

m Nk-1 yu 

Jn{Pn)< E ^^-^fe^l^M)- 

fc=0 £=0 
Nk>0 

This inequality holds, since for a spherical cap C we have 

N-1 m / ^ 

- - A(C) = 5; E^^ U ^ M-)-A(C) 

ra=0 fc=0 ^=0 \ tceQfe I 

Affc>0 ^ 
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Thus we can use Theorem 1141 to obtain 

MP.) < E E ^^^'-^"^^ = 4^ E E V = 

fc=o e=o k=o e=o ^k=o ^ 

Nk>0 Nk>0 

< 4^(. - 1) E ^ < 472^.-'"/= < 4V^^'7-1) ' . 

The estimate follows from the identity (a — l)(a + 1) = — 1. □ 
Corollary 16. 

(1) Let P/v be a (0, m, 2)-nei in base b and let = ^{Pn) ^ S^- Then the spherical cap 
discrepancy D{Zn) is bounded by 

D{Zn) < 44^/26" L™/2J. 

(2) Let P/v be the first N points of a (0, 2)-sequence in base b and let Z^ = $(P/v) C S^. 
Then the spherical cap discrepancy D{Zi\i) is bounded by 

D{Zn) < 44^2 (6^ + b^/A for all N. 



Note that Item (2) improves upon Theorem 11 by a factor of \/IogTogiV and hence, asymp- 
totically, spherical digital sequences are better than random sequences almost surely. 



m 


6 


7 


8 


9 


10 


11 


12 


13 


N = T^ 


64 


128 


256 


512 


1024 


2048 


4096 


8192 


log Af 


0.8829 
0.4329 


0.8436 
0.3829 


0.8279 
0.3515 


0.8632 
0.3456 


0.8518 
0.3235 


1.2128 
0.4392 


1.2285 
0.4259 


0.9546 
0.3180 


m 


14 


15 


16 


17 


18 


19 


20 


21 


N = 2"^ 


16384 


32768 


65536 


131072 


262144 


524288 


1048576 


2097152 


D(Zj^)*N^/^ 

VlogTV 
D{Z!^)*N^/^ 
log AT 


0.7925 
0.2544 


0.8862 
0.2748 


1.0331 
0.3102 


0.8337 
0.2428 


0.8562 
0.2424 


0.9854 
0.2715 


1.1167 
0.2999 


1.1463 
0.3004 



Table 1. The value D{Z]sr) denotes the maximum of the absolute values of 
the local discrepancies of Z]\f at the spherical caps centered at the spherical 
digital net Z^ based on a two-dimensional Sobol' point set. 



The numerical experiments shown in Table [T] seem to suggest that the correct order of the 
spherical cap discrepancy of spherical digital nets is 

— tor some 12 < c < 1. 

iV3/4 / - - 

In those experiments we calculated the spherical cap discrepancy for all spherical caps centered 
at a spherical digital net based on a two-dimensional Sobol' point set (this of course gives 
only a lower bound). 
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5.2. Fibonacci lattices. The Fibonacci numbers Ffn are given by -Fi = 1,^2 = 1 and 
Ffn = Fm-i + Fm-2 for all m > 2. A Fibonacci lattice is a point set of F^ points in [0, 1)^ 
given by 



\ Fjyi Fjyi J J 

where {x} = x — \ x\ denotes the fractional part for nonnegative real numbers x. The set 

'■= {fo, • • • , fpm-l} 

is called a Fibonacci lattice point set. 

The spherical Fibonacci lattice points are then given by 



m 1 



Zn = *(/„), < n < F„ 
and the point set 

Zf^ = {zq, . . .,Zf^-i} 

is the spherical Fibonacci lattice point set. 

In the following we prove a bound on the isotropic discrepancy of Fibonacci lattices, see 
also [351 ESI EZ]. 

Lemma 17. For the isotropic discrepancy Jp^ of a Fibonacci lattice J-m we have 



A:^2lFm if rn is odd, 

^^W\n if-' 



m IS even. 



Proof. Consider the case of odd integers m first. From [43\ Theorem 3] it follows that for 
m G N the Fibonacci lattice J-2m+i can be generated by the vectors 

0-2m+l = (-?^m/-F2m+l, ( — 1)™ ^ Fm+l / F2m+l) &2m+l = (-^m+l /-f2m+l , ( — 1)™'-Pm/-F2m+1 ) • 

This means that 

F2m+l = {ua2m+l + ?;&2m+i -.u^ve Z}n [0, 1)^ 

Let 

U{y) = {y + x e [0,1]2 : x = sa2m+i + th2m+i,^ < s,t < 1}. 

We call U{f^) a unit cell (belonging to the point Note that the area of a unit cell is 

^/F2m+i and each unit cell contains exactly one point of the lattice, see [33j . 

Since a2m+i -L ^'2m+i it follows that the minimum distance between points of the Fibonacci 
lattice is 



dmm{F2m+l) — niin{ || a2m+l || , ||&2m+l| 



F"^ ^ F"^ 

^ m ^ m+1 i 



F2m+1 V -^2m+l 

Thus the diameter of a unit cell is ■sJljF2m+\- 

Let now A be an arbitrary convex subset of [0, Let W° denote the union of all unit cells 
fully contained in A and let W denote the union of all unit cells with nonempty intersection 
with A or its boundary. 

We have 

U{fJ-X{A)<- lwifn)-HW) + X{W\A) = X{W\A) 



F2m+1 „ F2m+1 „ 

71=0 n=0 
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and 



1 

F2iri+1 



-F2m + 1 — 1 



F2r, 



U(/„)-A(A)>-^ ^ l^.(/j_A(t^°)-A(^\W^°) = -A(^\W^°). 

n=0 n=0 

Since the set A is convex, the length of the boundary of A is at most the circumference of 
the unit square, which is 4. Further we have 

W\A<Z{x(^[Q,lf\A: \\x - y\\ < y^2/F2m+i for some y e A} 

and therefore 

X{W \ A) < X{{x G [0, 1]^\A: \\x - y\\ < y^2/F2m+i for some y e A}) < 

where the last inequality follows since the outer boundary of enclosing set has length at most 
4 (which is the circumference of the square [0, 1]^). 
On the other hand 

A\W° Q {x e A : \\x - y\\ < y^2/F2m+i for some y G [0, 1]^ \ A} 

and therefore 



X{A \ W°) < X{{x eA:\\x- y\\ < for some y G [0, l]^ \ A}) < 4^2/^2^+1, 

since, by the convexity of A, the boundary of A has length at most 4 (which is the circum- 
ference of the square [0, 1]^). 
Thus we obtain 

-F2m + 1 — 1 



1 



2m+l 



U(/J-A(^) 



n=0 



< 4, 



F 



2m+l 



Now we consider even integers n = 2m with m > 2. Using the identity FmF2m-i 
Fm-iF2m = {-l)'^~^Fm, we obtain FmF2m-i = {-l)"'~^Fm (mod F2m)- Consequently 



k 



F2r, 



kF: 



2m- 1 



F2r 



F2m 



(-ir-^F„ 

F2m 



for k = Fm.. 



Thus the Fibonacci lattice has the equivalent generating vector 

0'2m 



, F2m F2m 

Analogously, using the equality Fm+iF2m-i- FmF2m = 
vector 

Fm+l ( — 1)™'-Pm-1 
F2m. ' F2m 

The area of the parallelogram spanned by a2m and 62m is 



l)™-Pm-i) we obtain the generating 



■'2m 



det 



FmFm—l ~\~ -^m-^'m+l 



1 



7?2 

-^2m 



F, 



2m 



Fm/F2m ( — I)'" Fm/F2m 
Fm+l/F2m {—^)^Fm-l/F2m 

Thus the parallelogram spanned by a2m and b2m does not contain any point of the Fibonacci 
lattice in its interior (i.e. is a unit cell of the Fibonacci lattice, see |15ii43j). Thus we have 



Let 



J^2m = {ua2m + vb2m : li, G Z} Pi [0, 1)^. 
U{y) = {y + x e [0,1]^ : x = sa2m+i + tb2m+i,0 < s,t < 1}. 
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We call U{f^) a unit cell (belonging to the point Note that the area of a unit cell is 

l/-^2m+i and each unit cell contains exactly one point of the lattice. 
Now 

2Fl 



1^ l|2 
|a2m|| - ^ 



and 



1712 

^2m 



lU l|2 _ +-^m-l ^ + Fm{2Fjn-l) 2F^ _ 2 
\02m\\ - ^ > ^ > - \\a2m\\ ■ 



2m 2m 2m 

Further it can be checked that \\a2m + ^^2™!!) ||a'2m ~ &2m|| > ||ct2m||- Thus the minimum 
distance between points of the Fibonacci lattice is 

, , ,, ,, V2Frn 
(J'm.m\-r2m) — l|<*2m|| — • 

^2m 

The diameter of a unit cell is given by ||a2m+''2m||- Using the relations F^ = Fm-i+Fjn-2, 
+ Pm-1 = F2m-i and F2m = (2-Pm-i + Fm)Fm we obtain 

^2m\\^2m + &2m|P = {Fm + -?^m+l)^ + {Fm — -^m-l)^ = 2F2m + 4-F2m-l + < 8F2m- 

Thus the diameter of a unit cell is bounded by 

\\o.2m + b2m\\ < 



F2m 

The result now follows by using the same arguments as in the previous case. □ 

Corollary 18. Let J>„ be a Fibonacci lattice and let Zp^ = ^{T^) C S^. Then the spherical 
cap discrepancy D{Zp^) is bounded by 



AAy/2/Fm ifm is odd, 
A'i^J'&/Fjn ifm is even. 



m 


17 


18 


19 


20 


21 


22 


23 


24 


^ m 


1597 


2584 


4181 


6765 


10946 


17711 


28657 


46368 


-v/log Fm 
b{Zp^)*FU^ 


0.6729 


0.6373 


0.6228 


0.6661 


0.6953 


0.6890 


0.7427 


0.6900 


0.2477 


0.2273 


0.2156 


0.2243 


0.2279 


0.2203 


0.2318 


0.2105 


m 


25 


26 


27 


28 


29 


30 


31 


32 


^ m 


75025 


121393 


196418 


317811 


514229 


832040 


1346269 


2178309 


b{Zp^)*F^^ 
Vlog Fm, 

b{ZFm)*FT 


0.6957 
0.2076 


0.7249 
0.2118 


0.7531 
0.2157 


0.7205 
0.2024 


0.8562 
0.2361 


0.7455 
0.2019 


0.7862 
0.2092 


0.8082 
0.2115 



Table 2. The value D{Zp^) denotes the maximum of the absolute values of 
the local discrepancies of Zp^, at the spherical caps centered at the spherical 
Fibonacci points Zp^ ■ 
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The numerical experiments shown in Table [2] seem to suggest that the correct order of the 
spherical cap discrepancy of spherical Fibonacci lattice points is of order 

(log Fm 



7^3/4 



for some 1/2 < c < 1. 



In those experiments we calculated the spherical cap discrepancy for all spherical caps centered 
at the spherical Fibonacci points (this of course gives only a lower bound) . 

6. Level curves of the distance function and their properties 

The Euclidean distance of two points on S^, given hy w = $(n, and z = $(a,T), can 
be written as 



zf = 2{l - w ■ z) = 2\l - {I - 2v) {I - 2t) -A^/{1-v)v{1-t)t cos{2tt{u - a)) 



2{v+t-2vt) 

The well-defined boundary curve of a spherical cap C{w,t) has the implicit representation 
(8) F{a,T):= \\w — z\\'^ — 2 [1 — t) = 0, where z = ^{a,T) moves on dC{w,t). 



(Here, y^2(1 — t) is the distance from the centre w to a point on the boundary of C{w,t).) 
Relation ([s]) describes a level curve Cud of the distance function \\w — $(tt, t)|| (for w fixed) in 
the parameter space which is the unit square, cf. Figure [2] For further references we record 
that for each w there are, in general, two exceptional levels 



rl = 2{l-t^) = \\w-pf = 4v, 

= 2 (1 - 4) = 2 (1 + t^) = \\w + pf = 4{l-v), 

where the boundary of the spherical cap centred at w passes through the North Pole (p) and 
the South Pole (— p), respectively, which may coincide if w is on the equator. For these level 
curves the singular behaviour at the poles imposed by the parametrisation $ plays a role. 

Suppose u = 1/2. Because the sign of the difference u — a is absorbed by the cosine 
function in the distance function, a level curve (a level set) is symmetric with respect to the 
vertical line a = u. The shape of the curves (sets) do not change when the point w is rotated 
about the polar axis except a part moving outside the left side of the unit square enters 
at the right side ("wrap around"). This "modulo 1" behaviour complicates considerations 
regarding convexity of the pre-image of a spherical cap centred at w under the distance 
function. Similarly and in the same sense, the level curves (sets) are symmetric with respect 
to the vertical line a = u it 1/2 mod 1 which passes through the parameter point of the 
antipodal point —w = $(u it 1/2 mod 1,1 — v). When identifying the left and right side 
of the unit square [0,1]^, we get the "cylindrical view". Clearly, all level curves except the 
critical ones associated with the distance to one of the poles are closed on the open cylinder. 
Thus, the critical level curves separate the open cylinder into three parts corresponding to 
the cases when neither pole is contained in the spherical cap centred at w, only one pole is 
contained in the spherical cap, and both poles are contained in the spherical cap. In both 
the first and last case the level curve can not escape the level set bounded by a critical curve 
(and the boundary of the cylinder). It is closed even in the unit square provided the level 
set is contained in its interior. In the middle case the level curves wrap around the cylinder; 
that is, start at the left side of the unit square and end at its right side at the same height; 
cf. Figure [2j 
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Let w be the North or the South Pole. Then the level curves are horizontal lines in 
the unit square which are smooth curves. Moreover, the pre-image of a spherical cap centred 
at one of these poles is a convex set (a rectangle). 

Let w be different from either pole (that is < v < 1). We use the signed curvature 
of the implicitly given level curve to determine the segments where it is convex (concave). 
First, we collect the partial derivatives up to second order: 

(9a) Fa{a, r) = — IGvr-y/ {1 — v)v {1 — t)t sin(27r(n — a)), 

(9b) Fria,T) = 4 ( l-2v-{l- 2r) / j^^ ~ ^ ll cos(2^(n - a)) J , 




(9c) Faa{a, t) = 327rV(l -v)v{1-t)t cos(27r(ti - a)). 



(9d) Faria, r) = -Svr (1 - 2r) " ^ sin(27r(u - a)), 



(9e) Frria, t) = -. ^ = cos(27r('u — aS). 

^ ' ^ ^ (l-r)r^(l-r)r ^ ^ " 

We observe that the partial derivatives involving differentiation with respect to r become 
singular as r approaches (North Pole) or 1 (South Pole). 
The signed curvature at a point (a, r) of Cw is given by 



(10) K = K(a, r) 



7? 7?2_27? /?/?_|_7? 7?2 



{F^ + F^f 



First, we discuss the denominator. Substituting the relations (9a) and (9b) we obtain 

F^ + F^ = 16 (l-2v-{l - 2r) cos(27r(n - a))] 

V V(l-r)T '7 

+ 16^2 {1-v)v{1-t)t [sin(27r(u - a))f . 

A necessary condition for the vanishing of F^ + is sin(27r(n — a)) = 0; that is, either 
a = uoTa = uzizl/2 mod 1. In the first case one has cos(27r(ii — a)) = 1 and, therefore, 
F^ + F^ = if and only if r = In the second case one has cos(27r(n — a)) = —1 and, 
therefore, F^ + F^ = if and only if r = 1 — It follows that the denominator of the 
curvature formula (10) vanishes if and only if z = $(a,r) (on the boundary of the spherical 
cap centred at w = ^{u,v)) coincides with w (that is the spherical cap degenerates to the 
point w) or z coincides with the antipodal point of w (that is the closed spherical cap is the 
whole sphere). In either of these cases the pre-image of the spherical cap under $ is convex. 

Suppose that the spherical cap is neither a point nor the whole sphere. Then 
F^ + F^ ^ 0. For the numerator in (10) we obtain after some simplifications 

FaaFj- 2FctTFaFr -\- F-j-q-F^ 

5127rV(l -v)v{1-t)t ( 1 - 2t; - (1 - 2r) ^^''^ ~ '^^ cos(27r(ti - a))] cos(27r(u - a)) 

v(i-^)^ / 




+ 5127r (1 — v) V — = cos(27r(n — q)) [sin(27r(n — a))Y 
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- 10247r2 (^1-v)v{1-2t) ll-2v- {1- 2r) ^''^ ~ ^ cos(27r(u - a))] [sin(27r(u - a))f 

V v(i-^)^ / 

The right-hand side above can be written as a polynomial in a; = cos(27r(n — a)) as follows 

FaaF^ - 2F^rFaFr + FrrF^ = Ax {1 - 2v - BHxf + ^ ' 



5127r2 (1 -r)r 
2ABH {l-2v- BHx) (l - x^) , 



where 



A = 5127rV(l -v)v{1-t) t, B = =, F = 1 - 2r. 

v(i-^)^ 



Reordering with respect to falling powers of x, we observe that the coefficient of vanishes 
and after simplifications we arrive at 



F,„F2 - 2F^rFaFr + F^rF^ = - AB (bH^ + ^) 



X^ 



+ A ( 2B^H + (1 - 2vf + ^t^^VtT ^ ) ^ - 2^^^ (1 - 2^) • 

OlZTT^ (1 — Tj T/ 



The coefficient of does not vanish for < f < 1 and < r < 1. Hence, we divide and get 

x^ +px + q=:Q{T; x) = :Q{x), 



TP — 2F F F -\- F 



AB ( BH^ 



+ 



5127r2{l-T)r^ 

where (using the definitions of A, B, and H) 

{l-2vf + B^l + 2H^) 2{1-2v){1-2t) 
P = p{r) = DO,-, , rr2^ ' <} = Q{V- 



B'^{i + H^) ' ' B{i + m) ■ 

We observe that p(l — r) = p{t) andq{l — T) = — ^(r). Hence, Q{t;x) = —Q{1 — t;—x) for all 
x. In particular, if ^ is a zero of Q{t; •), then so is — ^ a zero of Q{1 — r; •) and vice versa. The 
monic polynomial Q of degree 3 with real coefficients has either one or three real solutions 
(counting multiplicity). With the help of Mathematica we find that the discriminant of the 
polynomial Q is positive, 

discr(g) = -V - 27g2 > 0; 
that is, the polynomial Q has three distinct real roots. 
For t; = 1/2 the polynomial Q reduces to 

N 3 , X. 1 + 2H^ 1 + 2(1 -2t)^ 

Q{x) = x +px, where p = - ^ ^ =- ^ ^ _ ^^^2 ■ 

The solutions ±1 (if r = 1/2) correspond to $(Q,r) = ztiu and can be discarded, since we 
assumed that the spherical cap is neither a point nor the whole sphere. The solution zero 
yields that cos(27r(u — a)) = 0, which in turn shows that the zeros of the curvature (10) form 
the vertical lines at a = Mibl/4 mod 1 \i v = 1/2. 
Let V 7^ 1/2. Suppose Q has a zero at ±1. Then 



= Q(1) = l+p + q = ± 



B'^{i + m) B^{i + m) 
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Figure 2. Contour plot of the distance function \\w — $(a,T) || for various 
w. Level curves for distances ^/2 and 2^/v, 2y/l — v are emphasised. The solid 
curve shows vanishing curvature. 



which can only happen when t = v. This implies that $(a, r) = w, which is excluded by our 
assumptions. Suppose Q has a zero at 0. Since v / 1/2, this can only happen when r = 1/2. 

Having established that —1, (except when r = 1/2), and 1 cannot be zeros of the 
polynomial Q, we use Storm's theorem to show that the polynomial Q has precisely one 
solution either in the interval (—1, 0) or in the interval (0, 1) if r 7^ 1/2, cf. Tablejsj First, we 
generate the canonical Storm chain by applying Euclid's algorithm to Q and its derivative: 

Po(x) = Q{x) = + px + q, 

Pi{x) = Q'{x) = +p, 

P2ix) =pi{x)qoix) -po{x) = -yx-g, 

27g2 - V - 27q^ discr(Q) 
P3{x) = P2{x)qi{x) -pi{x) = -p- — = — = > 0, 

Pi{x) = 0. 

Let crlx) denote the number of sign changes (not counting a zero) in the sequence 

{po(x),pi{x),p2{x),ps{x)}. 

For X = we obtain the canonical Storm chain — (7, discr(Q)/(4p^)} and we conclude 
that (t(0) = 2 for (1/2 — v){l/2 — r) > and a{0) = 1 otherwise. For x = 1 we have 

(l-2v + BH)^ 
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Pl(l) =3+p 



2 + {l-2vf 



>0 iS {1 -t)t <3{1 -v)v, 



P2{1) 



2p 



B2 (1 + 2if2) + (1 - 2vy - 3BH (1 - 2v) 



3fi2 (1 + H^) 



> 0, 



f3(l) = ^^>0. 



(The positivity of P2(l) has been verified using Mathematica.) Hence, in ah three cases 
Pi{l) < 0, = 0, and pi{l) > one gets a{l) = 1. For x = — 1 we have 



^'o(-l) 



-l-p + ^ 



pi{-l) = 3+p 



(1 - 2V + BHY 
2 + H'^ (1 - 2t;)^ 



1 + i72 ^2 (1 + m) 



>0, 



>0 iff (1 - T)r < 3(1 - 



P2{-1) 



2p 
3 



^2 (1 + 2^2) + (1 - 2w)2 + 3BH (1 - 2v) 



3^2 (1 + H"^) 



<0, 



f3(-l)=^^>0^ 



Here, we obtain ct(— 1) = 2. Thus, by Storms Theorem, the difference o"(— 1) — ct(0) gives the 
number of real zeros of Q in the interval (—1,0] and a"(0) — is the number of zeros in 
(0,1], see Table |3} 

Table 3. The number of real zeros of Q in the intervals (—1, 0) and (0, 1) as 
it follows from Sturm's Theorem. 



Range of 


V and r 


a{-l) 


a{0) 


a{l) 


(-1,0) 


(0,1) 


Q<v < 1/2, 


< r < 1/2 


2 


2 


1 





1 


^<v < 1/2, 


1/2 < r < 1 


2 


1 


1 


1 





1/2 < v < 1, 


< u < 1/2 


2 


1 


1 


1 





1/2 < V < 1, 


1/2 < r < 1 


2 


2 


1 





1 



Having established that the polynomial Q has to every < r < 1 precisely one zero in the 
interval (—1,1) (cf. Table [s] and previous considerations), it follows that to each such zero 
X = x{t) there correspond two values of a by means of the trigonometric equation 

(11) cos(27r(n — a)) = a;(r). 

Because of the continuity of the coefficients in the polynomial Q (if < r < 1) the zero x(t) 
is also changing continuously and so are the solutions ai and 02- A jump can happen when 
they are taken modulo 1. We further record that along the vertical lines a = u it 1/4 mod 1 
one has 

2 (l-z;)t;(l-2T;)(l-2r) 



k{u±1/4,t) = -levr" 



[I- 4(1- v)v{l- 47r2 (1 - r) r)] 



3/2 
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which vanishes at r = 1/2 and along the vertical lines a = nib 1/2 mod 1 and a = u one has 

9 \/ (I — v) V (1 — t) T 

k{u, t) = n{u ± 1/2, r) = -Svr^ ^ ^ 

(1 - 2v) V(l - r) r + (1 - 2t) ^(1 -v)v 

which vanishes only as r — )• or r — >• 1 (if < ?; < 1). When identifying the left and right 



side of the unit square, these two lines separate the two solutions of ( 11 ) in such a way that in 
each part the points at which K,{a, r) vanishes form a connected curve varying about the "base 
lines" a = u lb 1/4 mod 1. It follows that these curves (together with the boundary of the 
cylinder) divide the cylinder into two parts in each of which the curvature K{a, r) has the same 
sign. The shapes of these curves do not change when w is rotated about the polar axis. We 
may fix u = 1/2 and because of the symmetries (including relation Q{t; x) = —Q{1 — r; — x)) 
it suffices to consider the curve of the zeros of k(q, r) for < v < 1/2 (recall that these curves 
are vertical lines for v = 1/2) and < r < 1/2 which lies in the strip < a < 1/2. We know 
that the zero x = x(r) of Q (we are interested in) in the given setting is in (0, 1) (cf. Table [3|. 
Using X to denote the derivative of x with respect to r, implicit differentiation gives 

(12) Q'{x{t)) x(t) = -p{t) x{t) - q{T), 

where it can be easily seen that Q'{x{t)) < 0, since x{t) is simple and Q{x) has a negative 
global minimum for positive x. For r's in (0, 1/2) at which x{t) vanishes one has 



(13) x{t) = -(?(r)/p(r) 

which follows by substituting 



l-6(l-r)rl-2v {l-v)v 
1 - 6 (1 - v) t; 1 - 2r V (1 - r) r ' 



= _ (l-6(l-.).)(l-2r) ^ (l-6(l-r)r)(l-2.) (l-r)r^ 

2{l-v)v{l-2{l-T)Tf 2 (1 - r) T (1 - 2 (1 - r) rf ^ {l-v)v 

For the second derivative of x at such r's we get 

I — 2v 

Q'{x{t)) x{t) = -p{t) x{t) - q{T) = -j^ — , . 

4 ((1 - r) rf/^ V(l -v)v[1-2t{2-t(?,- 2t))] 

The square-bracketed expression is strictly monotonically decreasing on (0, 1/2) and evaluates 
to zero at r = 1/2. Thus, the left-hand side has to be positive for all critical r in (0, 1/2) 
which in turn implies that x(r) < at such r's. We conclude that x{t) has a single maximum 



in (0, 1/2), since it cannot be constant as can be seen from (13) by evaluating x{t) at that r' 



at which 1 - 6 (1 - r') r' = and x{t) ^ oo as r ^ 1/2. By ([11| (recall u = 1/2) 

— cos(27ra(r)) = x(r) 
and it follows that a(r) has a single maximum in (0, 1/2). 

Proposition 19. The set of zeros of K{a, r) and the horizontal sides of the unit square divide 
the unit square either into three parts of equal sign of K{a,T) or in four parts, cf. Figure^ 

Using the trigonometric method, we obtain an explicit expression of the zero of Q in (—1, 1). 
The change of variable x = 2 y^—p/3 cos 9 gives the equivalence 

Q(x) = if and only if cos(30) = —\ — 

2p V —p 
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and therefore 



2V-p(r)/3 



cos - arccos , , , 
3 \2p{T) 





-1) 


l-pir)) 





(The discarded solutions are either smaller or larger than the given one. By our reasoning, 
they have to lie outside the interval (—1,1).) A limit process shows that x(r) — )• as r — )• 
or T — )• 1. Hence, when moving towards the upper or lower side of the square along the curve 
of zeros of the curvature, one approaches the corresponding "base line" a = ?x ± 1/4 mod 1. 
Eliminating the trigonometric term, along a level curve with parameter t (cf. ([8|) we have 



p p2 



-167r" 



(1 



1 1 



-8(l + 3ii^(l 



1 



r r 



(14) 



J? 7?2_2F F F -\- F F"^ 



2tY ] X + U{l-v)v{l- 2v) (1 - 2r) 



- 2G (3 + + 2,Gt - Zf) r 



where X = t— (1— 2t')(l— 2t). Reordering the terms and using the substitution G = t— (1— 2^;), 
we arrive at 

16 

'(l-r)V2 

+ 2 (2t + eG^t - 2t^ + 3G (3 - t^)) - 4 (5t - 3t^ + 3G (l + t^)) + 16tT^ 



The zeros of the numerator of the curvature ( |10[ ) are determined by a polynomial in r of 
degree 4. Thus, there can be at most four pairs (symmetry with respect to a = u) of points 
on the level set at which the curvature vanishes. For the sake of completeness, in a similar 
way one obtains 



(F2 + F2) 



■2n3/2 



1 



(G - 2tTf - 167r2 (1 - rf (G^ - 4Gtr - 4r (l - _ r)) 



3/2 



The critical curves when the boundary of the spherical cap passes through a pole are of 
particular interest. In the case of the North Pole (that is t = 1 — 2v, or equivalently G = 0), 
the curvature along the corresponding level curve reduces to 

,{1-2v){1-t)[2{1-v)v-{1 + 6{1-v)v)t + 2t^] 

(15) k{t] = —Ibir 



3/2 



(1 - 2vf - 167r2 (1 - r) V (r - 4 (1 - v) v) 

For given < r < 1, the corresponding value(s) of a can be obtained from the relation 

(1 - 2v)t 



(16) 



cos(27r(ii — a)) 



2^/(1 -v)v{1-t)t' 

From ([s]) it follows that the range for r is (0, ti] with ri = 4v{l — v). For future reference we 
record that for < v < 1 with v 1/2 the curvature ( |15| ) vanishes only for 

(17) : ^ 



(^l + 6{l-v)v- Vl - 4 (1 - d) (1 - 9 (1 - w) v)^ . 



(The other solution lies outside the interval [—1,1] as one can verify with Mathematica.) 



Proposition 20. The curves of zeros of the curvature (10) (as functions ofr) assume their 
extrema at r„ and 1 — Ty with Ty given in (17). 
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Proof. Suppose that u = 1/2 and < v < 1/2. Then < < 1/2. On observing that the 



right-hand side of (16) for r = is also the zero x{Ty) in the interval (0, 1) of the polynomial 



Q, it can be verified with the help of Mathematica that the right-hand side of (12) vanishes 
and therefore x{t^) = 0; that is, the zero x(r) is extremal at r = t„. Using the symmetry 
relation Q{t;x) = —Q{1 — t; —x), the zero x{t) is also extremal at r = 1 — Tt,. By means 



of ( |16| ) this translates into extrema of the curve of zeros of the curvature (10). A shift in u 
(rotation of w about the polar axis) does not change the shape of the level curves and the 
general result follows. □ 



Substituting (17) into the right-hand side of (16) gives the extremal value a zero of Q in 



(—1,1) can assume, also cf. (13): 



X[Ty) = , 



y 1 + 2 (1 - i;) + y 1 - (1 - i;) u (^9 (1 - 2vy - 
It can be shown that x{Ty) is a strictly monotonically decreasing function in v which is 



symmetric with respect to v = 1/2. Hence |a;(ri,)| < x{0'^) = l/\/2. Using this bound in (16) 
yields that |a — (ii±l/4)|<l/8 (when wrapping around) . 

When moving along the critical level curve towards the lower side of the unit square, which 
is associated with the North Pole, we have 

1 — 2v 

lim k(t) = -32tt^- p. (l-v)v, 0<v <l,v^ 1/2. 

T-s>o |1 - 2v\^ 

We conclude that the critical level curve with t = l — 2v (associated with the North Pole) has 
precisely one symmetric (in the cylindrical view) pair of intersection points with the curves 



of zeros of the curvature function (10) at r = r„ in the strip < r < 1. A similar result holds 



for the level curve associated with the South Pole {t = —(1 — 2v)). 

Proposition 21. Let < v < 1/2. Then the level curves with t in the range —(1 — 2v) < t < 
1 — 2v have precisely one symmetric (in the cylindrical view) pair of intersection points with 



the curve of zeros of the curvature function (10). For t in the ranges —1 < t < —(1 — 2v) or 



1 — 2v < t < 1 there are either no intersection points, one pair of tangential points, or two 
pairs. 

The analogue result holds for 1/2 < t < 1. (For f = 1/2 the level curves for distance \/2 
are the verticals at a = nibl/4 mod 1 and coincide with the curve of zeros of K{a, r) and 
also coincide with the critical curves. The other level curves have no intersections.) 

Proof. Without loss of generality assume that u = 1/2. We have already established that 
either critical level curve has precisely one pair of symmetric (with respect to a = u) inter- 
section point with the two curves of zero curvature about the base lines a = u± 1/4 at the 
values T = and r = 1 — r^. These parameter values also give the position of the extrema 
of the zero curves, cf. Figure [2} The left zero curve Z is increasing for r in (0, r^,), decreasing 
for r in (r^, 1 — t^) and increasing again for r in (1 — r^,, 1). 

Let —(1 — 2v) < t < 1 — 2v and F^ denote the left half of the corresponding level curve 
starting at the left side at some point (0, ri) and ending at some point {u, T2). (The other half 
is symmetric.) We note that the part where the zero curve is increasing is contained in the 
regions separated off by the critical level curves. Thus, an intersection between Z and F^ can 
only occur for r in the interval [ty,l — ty]. The curvature along F^ changes continuously (cf. 
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(14) and subsequent formula) from negative to positive value. Hence, there is an intersection 
point of Z and and the Yt cannot change abruptly. In particular, both Z and Fc]^ pass 
through (1/4, 1/2), which is their only intersection point because for t ^ 1/2 the vertical line 
a = 1/4 separates both curves. A Ft with 0<t< 1 — 2v (—(1 — 2v) < t < 0) has to intersect 
Z in the strip l/4<Q!<l/2(0<a< 1/4). Inspecting the partial derivatives of F (cf. 
([9])) it follows that the Gradient of F at the intersection point, which is the outward normal 
at the level curve Ff, points into the upper left part; that is the tangent vector at Ff at the 
intersection point shows to the right whereas the tangent vector at Z at this point shows to 
the left. Moreover, if < t < 1 — 2v, then the curve Fq separates Z and F^ for r > 1/2 and in 
the remaining part both curves Z and Ff bend away from each other because Z is decreasing 
with growing r and the curvature along Ff becomes positive. Consequently, there is only one 
intersection point of Ff and Z. A similar argument holds for —(1 — 2v) < t < 0. 

Let 1 — 2f < t < 1. Let Ff denote the left half of the level curve. For t sufficiently close 
to 1 there is no intersection with Z and the level curve is convex. If Ff and Z intersect in 
only one point, then the level curve is still convex, since the curvature function K{a, r) has 
positive sign in the section between Z and the vertical line a = u. In this case Z and Ff share 
a common tangent at the intersection point. If the curvature along Ff changes its sign to 
negative, then it has to become positive again, since it is positive when crossing the vertical 
a = u. But after changing back to positive curvature, both curves are bending away from 
each other. So, there can be no other intersection point. 

By symmetry with respect to the line a = u one has pairs of symmetric intersection points. 

Shifting u does not change the form of the curves and their relative positions. This com- 
pletes the proof. □ 

7. Proofs 

Proof of Lemma For t = 1 the spherical cap is a point and for f = — 1 it is the whole 
sphere. Their pre-images (a point and the whole unit square) are convex. So, we may assume 
that -1 < t < 1. 

Case (i): Let w be either the North or the South Pole. Then the pre-images of the 
boundary of spherical caps centred at w are horizontal lines in the unit square. Hence, the 
pre-image of such a spherical cap is convex. 



Case (ii): Let w be on the equator (that is v = 1/2). We know that the curvature (10) 
vanishes along the lines a = u ±1/4 mod 1. First, suppose that u = 1/4. Then the pre- 
image of any spherical cap centred at w with boundary points at most Euclidean distance 
V2 away from w is convex. For a larger spherical cap C it follows that its complement C 
with respect to the sphere (centred at the antipodal point —w) has the property that points 
on the boundary have distance < \/2 from —w. Hence, the pre-image of C is convex. When 
rotating w about the polar axis (that is shifting n), the vertical boundaries of the square cut 
these convex sets into two parts. We conclude that the pre-image of a spherical cap centred 
at w or its complement with respect to the sphere is the union of at most two convex sets. 

Case (iii): Let w be neither the poles nor located at the equator. Without loss of generality 
we may assume that w is in the upper half of the sphere; that is < u < 1/2. (Otherwise we 
can use reflection with respect to the equator.) First, let us consider the canonical position 



1/2. Let —(1 — 2v) < t < 1 — 2v. Then, by Proposition 21, there are precisely two 



(symmetric) points along the level curve at which the curvature vanishes, say at (ai,Tf) and 



^The spherical cap is the half-sphere. 
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{02, Tt). This yields a decomposition of the unit square into three vertical rectangles such 
that either the part above or below the level curve is convex. Let 1 — 2v < t < 1. By 
Proposition 21 the level curve is already convex or there are two pair of symmetric points at 
which the curvature along the level curve vanishes and a sign change occurs. Hence, there 
are numbers ri < T2 such that the level curve is convex for r < ti and convex for t > t2- 
The remaining middle part can be covered by a convex isosceles trapezoid which in turn can 
be split by some vertical line contained in the level set associated with the level curve. Thus, 
one has again two convex polygons which are divided into a convex and non-convex part by 
the level curve. A similar argument holds for —1 < t < —(1 — 2v). 

A shift of u does not increase the number of vertical rectangles needed for < t < 1 — 2f . (In 
fact, one may even reduce the number of elements of the partition.) In the case 1 — 2t; < t < 1 
one may need to use a covering of the pre-image of the spherical cap with up to 7 pieces. A 
more precise analysis is listed in Table |4j □ 



Table 4. Worst-case admissible convex covering with p part and q of which 
are convex. The vertical lines show canonical positions of the vertical borders 

of [0,1]2. 





p q 


2p-q 


A 


4 


2 


6 


B 


5 


3 


7 


C 


7 


3 


11 • 


D 


6 


3 


9 


E 


6 


4 


8 




Proof of Proposition^^ Radon's theorem (see e.g. |5l Theorem 4.1]) states that any set of 
d + 2 points from M*^ can be partitioned into two disjoint subsets whose convex hulls intersect. 
Particularly, let A denote a set of 5 points on the sphere. Then by Radon's theorem there 
exists a partitioning of A into disjoint subsets whose convex hulls intersect. Thus the set A 
can not be shattered by the class of halfspaces. Since every spherical cap is the intersection 
of the sphere with an appropriate halfspace, the set A can also not be shattered by the class 
of spherical caps. Thus the VC dimension of the class of spherical caps is at most 5. 
On the other hand, let the set A consist of the points of a regular simplex, which lie on the 
sphere. Then some simple considerations show that the set A is shattered by the class of 
spherical caps. Thus the VC dimension of the class of spherical caps (and therefore of course 
also the VC dimension of the class C of spherical caps for which the centre w and the height 
t are rational numbers, which was used in Section |4]) equals 5. □ 

Proof of Theorem As mentioned directly after the statement of Theorem[9j the lower bound 
in the theorem follows directly from ^ . To prove the upper bound we use Theorem [Tj By 
Proposition [S] the VC dimension of the class C in Section [4] is 5. For simplicity we assume 
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that the constant in Theorem [7| is an integer. Then for any s > 2K 



D{Zn) > 



< 



1 fKs"' 



2\ 5 



and consequently we have 

E(D(Z^)) < 



2K 



{s+l) 



< 



< 



s=2K 
oo 

+ E . 



N 



D{Zn) > 



2K 



s + l fKs^\^ _2,2 



K 



N 



for some appropriate constant K. This proves the theorem. 



□ 



Proof of Theorem 10. Let C* C C denote a hemisphere, i.e. a spherical cap whose normalised 
surface area measure is c(C*) = 1/2. By the central limit theorem for any t > 



N-l 



n=0 

and consequently, for any given e > and sufficiently small C^{e) > 0, 



du as N ^ oo, 



f 1 

I n=0 



for sufficiently large N. Since D{Zn) > jr lc*{Xn) - ^(C*) 



(18) 



'{d{Zn)<C3N-'/^] Kef. 



this implies 



for sufficiently large A^. 

On the other hand, by Theorem [7] for any given e > and sufficiently large C4{e) 



(19) 



D{Zn) > 



< e/2 



for sufficiently large A^. Combining (18) and (19) we obtain 

C3 



< D{Zn) < 



> 1-e 



for sufficiently large N, which proves Theorem 11 



□ 
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